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Abstract 

Many sources of errors exist in the manufacturing process of complex shapes. Some approximations occur at each step 
from the design geometry to the machined part. 

The aim of the paper is to present a method to evaluate the effect of high speed and high dynamic load on volumetric 
errors at the tool center point. 

The interpolator output signals and the machine encoder signals are recorded and compared to evaluate the contouring 
errors resulting from each axis follow-up error. The machine encoder signals are also compared to the actual tool center 
point position as recorded with a non-contact measuring instrument called CapBall to evaluate the total geometric errors. 
The novelty of the work lies in the method that is proposed to decompose the geometric errors in two categories: the 
quasi-static geometric errors independent from the speed of the trajectory and the dynamic geometric errors, dependent 
i on the programmed feed rate and resulting from the machine structure deflection during the acceleration of its axes. 

The evolution of the respective contributions for contouring errors, quasi-static geometric errors and dynamic geomet- 
ric errors is experimentally evaluated and a relation between programmed feed rate and dynamic errors is highlighted. 

Keywords: Machine tool, geometric errors, servo errors, dynamic errors, high-speed machining 
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1. Introduction 

High speed machining is known to reduce machining 
time and improve surface quality, due to its particular cut- 
ting process pQ . The ability to control the tool orientation 
with respect to the workpiece also leads to productivity 
improvements. 

The manufacturing process includes several steps. First, 
the part to machine, generally described by a digital model, 
is used to generate the appropriate tool path for machin- 
ing via a computer-aided machining software. This step 
generally comes with approximations in the part geome- 
try: a complex shape can be interpolated and discretized, 
leading to differences between the generated tool path and 
the numerical model of the part. 

The post processing of the computer-aided machining 
file is required to convert the tool path into CNC com- 
mands. The tool path, originally described in the work- 
piece coordinate system (WCS) must be transformed in 
the joint coordinate system (JCS). This inverse kinematic 
transformation is most often performed during the post 
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processing or sometime directly calculated by the NC-unit, 
but in any case, approximations may occur [TJ|2]. 

Then, the NC-unit interpolates the given tool path ex- 
pressed in the JCS to provide synchronous and dynami- 
cally admissible input for the machine drives. Each ma- 
chine axis is subject to follow-up errors, leading to orien- 
tation and position errors of the tool pQ . 

The machine structure, in its quasi-static state, suffers 
from link and motion errors causing position and orienta- 
tion errors of the tool [3-5J. Moreover, thermal variations 
of the machine is another source of errors [6] . 

When programmed feed rate and actual velocity in- 
crease, as is usual during the high-speed machining pro- 
cess, the dynamic solicitations of the structure become 
higher because of the acceleration of the different parts 
of the machine. This may result in some alterations of 
the machine geometry, causing further tool position and 
orientation errors pj [5] . 

Finally, the forces appearing during the cutting process 
lead to tool deflection, causing another source of machined 
surface location error [6]. 

Schmitz et al. compared the contributions of the geo- 
metric, thermal, controller and cutting forces errors to in- 
accuracies on a three-axis geometry, using different meth- 
ods and instruments to measure and evaluate the respec- 
tive contributions at different programmed feed rates |S]. 
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Zhu et al. identified the repeatable errors of machine tool 
using B-spline models using a cross-grid scale system [5]. 
Ibaraki et al. identified kinematic errors of a five axis ma- 
chine tool from geometric errors of finished workpieces and 
compared the results to ball bar measurements [10] . Bring- 
mann and Knapp measured the relative tool center point 
(TCP) deviations between a sphere and a nest of four lin- 
ear probes to calibrate a five axis machine tool jTTj . Erkan 
and Mayer used the machine own touch trigger probe to 
measure a set of sphere of uncalibrated relative position 
to calibrate the machine and assess its positioning per- 
formance [T3 . Finally, Schwcnkc et al. provide a compre- 
hensive overview of other measurement method to identify 
machine tool error in |14) . 

A method to evaluate the contributions of error sources 
to the Cartesian volumetric errors at the tool center point 
of a five axis machine tool at different programmed feed 
rates is described in this paper. The Cartesian volumetric 
errors are defined as the three components of the vector 
from the theoretical position of the TCP relative to the 
workpiece frame to its actual position. They are decom- 
posed into the contribution of three error sources: 

• the effect of follow-up errors of the axis drives, later 
called contouring errors and written S c ; 

• the quasi-static geometric errors of the machine, which 
include link and motion errors and thermal drift, 
written S qs ; 

• the dynamic geometric errors, resulting from the ma- 
chine structure deflection under dynamic load, and 
written d ( i- 

The evaluation of those three error sources can help 
quantifying the relative impact of dynamic errors to the 
total volumetric errors, and the relevance of an associated 
model for corrective actions. 

The paper begins with a presentation of the principle of 
the method is presented. Then, the mathematical models 
and the hypothesis drawn to evaluate each error source 
is developed. The experimental setup is described and 
finally, the results gathered from low to high speed motion 
are given and discussed. 



2. Principle of the method 

The nominal position of the tool center point relative 
to the workpiece is evaluated with the nominal direct kine- 
matic transformation (DKT) of the machine fed with the 
controller inputs of the machine. 

The actual position of the tool relative to the workpiece 
is measured with a non-contact sensor called CapBall, de- 
veloped in house [5] , which provides the position of a ball 
mounted on the machine table relative to the TCP. 

The difference between the nominal and the actual po- 
sition of the tool relative to the workpiece, pictured on 
Fig.[l] is due to the servo-errors of each axis of the machine 
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Fig. 1: Difference between the nominal position of the TCP relatively 
to workpiece x nom calculated from controller inputs and direct kine- 
matic model and its actual position t measured with the CapBall. 



and to all the geometric defects of the machine structure. 
This difference corresponds to the sum of S c , S qs and Sd, 
respectively the contouring error, the quasi-static geomet- 
ric error and the dynamic geometric error. 

The contouring error S c can be evaluated by recording 
the machine controller inputs and encoder actual valuesj 
and comparing the relative position of the tool and the 
workpiece from the two points of view. 

At low velocity, the dynamic geometric errors Sd can be 
considered negligible because the dynamic load is low with 
respect to the usual high stiffness of machines. Thus, once 
contouring errors have been accounted for, the remaining 
errors measured at a low programmed feed rate can be 
defined as the quasi-static geometric errors 8 qs . 

Finally when the programmed feed rate and the ma- 
chine velocity increase, alterations of the machine struc- 
ture under dynamic load cause the dynamic geometric er- 
rors Sd- They can be evaluated by measuring the total 
geometric errors along a single trajectory at different pro- 
grammed feed rates and removing the contouring error and 
quasi-static geometric error contributions. 



3. Model for error contributions 

3.1. Decomposition of the total errors 

During a measurement sequence, n points, equally dis- 
tributed, are recorded on each of the three channels of the 
CapBall. 



3 From position feedback linear encoders for linear axes and an- 
gular position feedback encoders for rotary axes. 
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Fig. 2: Decomposition of the total error into contouring error (subsection |3,2[ l, quasi-static geometric error (subsection |3-3[ > and dynamic 
geometric error (subsection 13. 4 J . 



The CapBall readings provide the values of the 3 com- 
ponents of the measured difference between P t - the actual 
tool center point - and P w - the actual center of the mas- 
ter ball mounted on the machine table -, written r for 
each of the n recorded points (see Fig. [T]): 



T = Pt-P w 



(1) 



The results are gathered in a n x 3 matrix written x, 
as shown in eq.pl): 



T\ 



(2) 



where tj, is the CapBall reading for fc-th point of the tra- 
jectory. The effect of the master ball and sensing head 
set-up positioning errors on volumetric errors are previ- 
ously identified with the method described in [12] and 
removed from \. 

The direct kinematic transformation (DKT) associated 
with the nominal machine geometry and the five axis con- 
troller inputs allows to calculate: 

• the position of Pt t nom ~ the nominal TCP - relative 
to the machine frame, considering the machine at its 
nominal geometry; 

• the position of P w _ nom - the nominal master ball 
center - also relative to the machine frame and also 
considering the machine with no geometric defect. 

The difference of position between those two points can 
be expressed as T nom calculated with the axis controller 
inputs (see Fig. [T] and eq.|3])). 



1-JJ\ J- (t> noun 7 +*-nomi -* nomi ^nom : ^n 



(4) 



The nominal difference matrix %,. 



calculated from 



the controller inputs - is formed as described in eq.([5]). 
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Then, the difference between x and Xnom expresses the 
effect of the contouring error, the quasi-static geometric 
error and the dynamic geometric error, as explicited in 
Fig. [2] and eq.^: 
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3.2. Contouring error 

The actual effect of follow-up errors on Cartesian vol- 
umetric errors is defined by the contouring error e, which 
is the orthogonal deviation of the actual tool path with 
respect to the controller inputs tool path in the part co- 
ordinate system |15j . The difference between Cartesian 
follow-up error and contouring error is schematically de- 
picted in Fig. [4} in case 1, the follow-up errors along x and 
y are smaller than in case 2, but the resulting contouring 
error ei is higher than 62- 

The difference T enc is calculated for each point, con- 
sidering the direct kinematic transform of the machine at 
its nominal geometry and the encoder values: 
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Fig. 3: Example of the three components of the Cartesian volumetric difference projected in the machine frame, computed with synchronised 
controllers inputs (xnom) and actual encoders values (xenc) for a programmed feed rate F = 5, 000 mm/rain. 
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Fig. 4: Schematic example of the difference between axis by axis 
Cartesian follow-up error 5 and contouring error e (adapted from 

ESI). 
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The encoder actual values volumetric difference matrix 
is built by the same process as in eq.pH): 



X.enc 




(9) 



The contouring error is evaluated by computing the dif- 
ference between the nominal volumetric difference matrix 
Xnom an d the encoder actual values volumetric difference 
matrix x enc - 

The controller input signals must be synchronised with 
the encoders actual values to evaluate the effect of the 
follow-up errors on the volumetric errors at the tool tip 
without taking into account the delay introduced by the 
recording process of the axis controllers signals. The record- 
ings performed show that the delay between encoder sig- 
nals and controller inputs signals was f 8 ms regardless of 
the programmed feed rate F. The encoder volumetric dif- 
ference matrix Xenc ^ s obtained by bringing backward the 
encoder values by 18 ms. The synchronised controller in- 
puts and encoders volumetric errors are plotted on Fig. [3j 

The contouring error along the experimental trajectory 



is given by eq.(10): 



(10) 




Fig. 5: Influence of the eight link errors on the geometric errors for 
the five axis machine tool studied. 



where S c has a n x 3 size, representing the contouring error 
expressed in the machine frame. 

Combining eqs.(|6]) and ( 10 ) allows to express the quasi- 
static geometric error and the dynamic geometric error as 



the difference of the recorded signals using eq.( 11 1 
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3.3. Quasi-static geometric errors 

3.3.1. Definition 

The quasi-static geometric error S qs is defined as the 
geometric error independent of the machine velocity. It is 
decomposed into three sources: 

• the link errors that represents the axis to axis loca- 
tion errors of the machine; 

• the motion errors of each axis; 

• the thermal drift of the machine. 

Consequently the quasi-static geometric error is writ- 
ten as: 

S qs = Si + S m + S td (12) 

where Si, S m and Std are respectively the contribution of 
the link errors, the motion errors and the thermal drift. 

3.3.2. Link errors 

Mir et al. proposed a model to express the effect of 
axis-to-axis link errors on volumetric errors [3]. In this 
model, the link errors for a five axis machine tool are de- 
scribed by a set of eight geometric errors of the machine 
structure. For the example in Fig. [5j the link errors are: 



• 5jYt the out-of-squareness between X and Y; 

• 5az, the out-of-squareness between Y and Z; 

• 50z, the out-of-squareness between X and Z; 

• 5 Pa, the tilt of A around y; 

• 5ja, the tilt of A around z; 

• Sac-, the tilt of C around x; 

• 5fic j the tilt of C around y; 

• dye, the offset of A relative to C in y. 
They are gathered in an array written Sq^. 

s Qi = (hv Sa z 5/3 z S/3 A 8j A 5a c Sf3 c Sy c ) 



(13) 



The model is based on a Jacobian principle. The vol- 
umetric errors due to the eight link errors of the machine 
at a pose k, written as a three components column vector 



Ti t k, is given by eq.(14| 
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where J^k is the link errors Jacobian matrix expressed at 
a pose k. More details are available in [HI [5]. Zargarbashi 
and Mayer described a method using CapBall data to iden- 
tify the eight machine errors 5, 12 . This method is also 
used here to identify the eight parameters of the model. 

Then, for each known pose of the machine, correspond- 
ing to the encoder actual values, the link errors Jacobian 
matrix Jjj, and then Ti f. are calculated, leading to the ex- 



pression for the contributions of the link errors Si (eq.( 15 )): 



S,= 
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Vl,n J 



(15) 



where Si has a n x 3 size, representing the three Cartesian 
components of the modelled contribution of the link errors 
to the volumetric errors at each point along the trajectory. 

3.3.3. Motion errors 

A reference thermal state of the machine, in which the 
value of the drift is considered to be zero, must be chosen 
for a test at low programmed feed rate. Then, the portion 
of measured errors during this test not explained by the 



link errors is calculated by cq.( 16 1 and is attributed to the 



motion errors and local dynamic errors. 
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(16) 



In [4], Bringmann and Knapp proposed a Fourier se- 
ries to model the motion errors of a link. Slamani et al. 
proposed a polynomial model for motion errors in |16j . 
In both cases the effects of motion errors are defined by 
continuous and smooth mathematical functions. 
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Fig. 6: Comparison of the local polynomial model for motion errors contributions <5 m and measured residual errors S m ,exp (sum of the motion 
errors and local dynamic errors) along a trajectory at a programmed feed rate of F = 1 000 mm/min. 



In order to evaluate the volumetric contribution of mo- 
tion errors, a polynomial model is used. The only aim of 
this model is to locally predict the effect of the five axes 
motion errors on volumetric errors along the experimen- 
tal trajectory to subtract it from the measured volumetric 
errors. The polynomial method proved to be easier to im- 
plement while providing satisfactory results. 

The dynamic geometric error proved to be negligible 
at a low programmed feed rate (compared to link and mo- 
tion errors) , and they are visible only localljP] along the 
trajectory, as depicted in Fig. [6] 

A normalised path parameter t n is created to describe 
the experimental trajectory. It has n linearly distributed 
components, from for the first point to 1 for the last 
point. The residual errors matrix is curve-fitted with poly- 
nomial functions of t n . The approximation leads to 3 iden- 
tified polynomials P x , P y and P z . The applied polyno- 
mial curve fitting leads to a filtering of local dynamic errors 
and measurement noise. 

Experimental results show that, for the machine tested, 
20 was a suitable value for the degree of the polynomial. 
Below 20, the curve fitting did not represent the residual 



4 Dynamic geometric error mainly appears at the beginning and 
the end of a trajectory, and at sharp corners. 



errors well, and above 20, some unexpected oscillations 
appeared. The model predicted contribution of motion 
errors is compared to actual measurement in Fig. [6] 

The contributions of motion errors to volumetric errors 



are modelled in eq.(17) 
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where t n ^ is the fc-th component of the normalised path 
parameter. 

The main weakness of this method is to compute mo- 
tion error contributions based on the time instead of the 
actual pose of the machine at each point. However, the 
small differential rate of the contributions of the motion 
errors reduces the effect of the approximation used. 

3.3.4. Thermal drift 

Even if a warm-up is performed and care is taken to 
favour the thermal stability of the machine, a small ther- 
mal drift can be present. 



As the execution of an experimental trajectory can take 
less than 15 s, it can be considered as instantaneous com- 
pared to the movement of the sensing head due to thermal 
drift. From that point of view, the effect of the thermal 
drift on the volumetric errors can be modelled as an offset 
from the reference thermal state on the measured volumet- 
ric errors. For each test, this offset is given by the mean 
value of the difference between the measured deviation x 
and the encoder only prediction with a nominal model, the 
modelled link errors and the motion errors contributions, 



as expressed in eq.(18) 
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(18) 



where % x is the i-th line of the matrix X- 

Finally, the drift, defined as constant for all the mea- 
surement point of a test is expressed by a n x 3 matrix 5td, 
composed of n lines equal to td: 



Std = 



f td\ 



,tdi 



(19) 



(rax 3) 



3.4- Dynamic geometric errors 

The dynamic geometric errors are here defined as the 
additional errors occurring when programmed feed rate, 
and so dynamic forces on the machine structure, increases. 
The dynamic errors result from varying alterations of the 
machine components under dynamic forces such as deflec- 
tion of the machine structure. 

Eqs.((6]) and (12) allow to express the dynamic geomet- 
ric errors Sd as the contribution that does not show at low 



speed, in eq.(20): 



Sd = (X ~ Xe 
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(20) 



4. Experimental setup and conditions 



4-1- Machine tool 

This study was carried out on a Huron KX8-five five 
axis machine tool depicted in Fig. [7| This machine has a 
WCAYFXZT structure, with a 45-degree-tilted A rotary 
axis. It is equipped with a Siemens Sinumcrik 840D Pow- 
er line numerical command unit. The NC unit allows the 
recording of system variables such as axis controller inputs 
and encoder actual values for each axis (X,Y,Z,A and C). 
The machine documentation specifies that both axis con- 
troller inputs and encoder actual values already include 
backlash compensations [17], so this potential source of 
error is not included in this study. 

Previous tests show that the axis controller inputs re- 
main the same for each execution of a program. Thus, the 
recording of controller inputs can be performed separately. 




Fig. 7: Huron KX8-five five axis machine tool fitted with the modified 
CapBall system. 



The encoder actual position and the CapBall measure- 
ments are recorded together during the execution of tra- 
jectory . 

Finally, all 10 channels (5 for the controller input and 
5 for the encoder actual values) are gathered. 

The NC cycle time is 3 ms and the sampling frequency 
of the CapBall is 10 kHz. To accurately synchronise and 
compare all the signals, the 10 NC signals are linearly in- 
terpolated at 10 kHz. A point is added at the beginning 
of the experimental trajectory to provide a tag dedicated 
to synchronisation of all signals together. More details are 
given at the end of subsection |4.3| 



4-2. External measuring instrument 

The CapBall consists of a precision master ball mounted 
on the table and a sensing head fitted with three capaci- 
tive sensors, mounted in a tool holder in the spindle (see 
Fig.0. 

The three capacitive sensors axes are nominally orthog- 
onal to each other and intersecting in one point called Pt, 
chosen as the virtual TCP. This forms an orthogonal mea- 
surement frame. 

The capacitive sensors are pre-calibrated for the spher- 
ical target using a set of three high precision linear stages. 
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Fig. 8: Experimental trajectory: left, evolution of each axis setpoint as a function of a normalised path parameter u representing the linearly 
scaled curvilinear abscissa and right, view of the master ball center P w trajectory in the machine coordinate system (attached to the Y-stagc). 



The response non-linearity is kept under 0.5% as long as 
the eccentricity is kept within ±300 \im and the distance 
between the sensor and the ball is in a ±300 fxm range 
centred on the position at which the sensor gives a Volt 
signal. The machine is programmed to keep the center 
of the master ball P w close to P< to avoid exceeding the 
measurement range of the sensors. In this range, stan- 
dard type A uncertainty on the sensors output has been 
evaluated below 0.4 \xm [T2] . 

The spindle orientation is regulated to keep the sensing 
head orientation constant relative to the machine coordi- 
nate system - or machine frame - (x, y, z) (Fig. [7]). 

To express the position of Pt relative to P w in the ma- 
chine tool frame, the orientation of the sensors axes are 
calibrated. Then, the position of Pt relatively to P w can 
be projected in the machine tool frame, to express t. 

The sensing head described in jS] was completely re- 
designed to significantly increase its stiffness and to ensure 
negligible displacement of the sensors during a measure- 
ment under dynamic solicitations in the maximum oper- 
ating range of the machine tool. The master ball base is 
a 40 mm diameter steel bar screwed to the machine table 
via a solid post (Fig. [7]). 

The sensors readings are acquired trough a LABView 
application. All the computation are implemented in MAT- 
LAB programs. 

4-3. Experimental trajectory 

To evaluate the contributions of the contouring, quasi- 
static and dynamic errors at different feed rates, an exper- 



imental trajectory is chosen to provide several character- 
istics. 

It must not last longer than 15 seconds when executed 
at the lowest programmed feed rate (i.e. lOOOmm/min), 
in order to avoid exceeding the NC unit recording buffer 
size. This constraint reduces the maximum length of the 
trajectory, forbidding performing standard tests as Schmitz 
et al. did in [6]. 

The motion of the master ball results from the rotary 
axis motion so its path is a complex curve. As the sensing 
head motion was programmed with linear interpolation, a 
resulting programmed difference of position between P w 
and Pt appears. The maximum admissible difference was 
chosen as 100 ^m in each direction of the machine frame. 

The experimental trajectory - with a nearly square cor- 
ner - was chosen to provide sizeable geometric discontinu- 
ities which require accelerations of several axes. This is 
the solution chosen to generate a large range of dynamic 
loads while the programmed feed rate varies. The cho- 
sen experimental trajectory is described for each axis as 
a function of a normalised parameter u varying from to 
1 in Fig. [8] This parameter represents the linearly scaled 
curvilinear abscissa along the trajectory. The trajectory, 
as perceived in the coordinate system attached to the Y- 
stage, is also depicted. The trajectory is described by 17 
linear segments (programmed with G01 blocks). The lin- 
ear interpolation format leads to acceleration peaks at the 



transition between each segment, as shown in Fig. 10 



The motion of each axis is programmed with linear in- 
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Fig. 9: Computed contributions of the follow-up errors, quasi-static and dynamic errors and the associated machine axes acceleration and 
block transition at a programmed feed rate of F = 1 000 mm/min. 



terpolation, and all the axes are synchronised (FGR0UP (X , 
Y, Z, A, C)). The machine is commanded in joint space 
(TRAF00F option)to provide the sought dynamic solicita- 
tions. To avoid critical speed decreasing at the corner 
formed by the nominal path at each block transition, a 
contouring option is used, allowing rounding the corners 
with an axial tolerance set at 10 fim. 

The trajectory is run at six different programmed feed 
rates F, geometrically distributed from 1 000 mm/min to 
18 896 mm/min. 

5. Results 

5.1. Mean value of the norm 

The mean percentage value of the norm along the ex- 
perimental trajectory, called 5k % for the source k, is used 



to compare the contributions of each source of error and 
is calculated as follows: 



>k,% 



E 



>qs\ 



***! 



x 100 



(21) 



where l 5u is the i-th line of the matrix dk, and k can 
be c for contouring errors, qs for quasi-static geometric 
errors and d for dynamic errors. Thus, 5k % gives the mean 
proportion of the total error due to the source k along the 
trajectory. The repartition of the error sources for the six 
tests is given in Table 1. 

For the geometric errors, the distinction is made be- 
tween the three identified sources. The experiments were 
carried out with deactivated compensation tables, avoid- 



0.5 0.6 

Time (s) 



0.06 
0.04 



0.02 



0.1 



x 10" 



0.2 



0.3 



0.4 



0.5 0.6 

Time (s) 



0.7 



0.5 0.6 

Time (s) 



0.8 



0.9 




0.8 0.9 




0.5 0.( 

Time (s) 



1 
0.5 



0.1 



0.2 



0.3 



0.4 



0.5 0.6 

Time (s) 



0.7 



0.8 



0.9 



Fig. 10: Computed contributions of the follow-up errors, quasi-static and dynamic errors and the associated machine axes acceleration and 
block transition at a programmed feed rate of F = 18 896 mm/min. 



ing any effect of software compensation of the link and 
motion errors. In an industrial context, the effect of link 
errors can be compensated, as it is a model-based calcu- 
lated contribution. 

The increasing relative impact of the thermal drift er- 
rors (see Table 1) may not be related to the increasing 
feed rate, but the drift is slowly increasing from test to 
test while the thermal condition are getting away of the 
reference thermal state. Moreover, the effect of the ther- 
mal drift observed during this series of tests is not reflect- 
ing actual machining conditions because the spindle is not 
running but only regulating its angular position. More 
sizeable errors can be expected under heavier spindle use. 

Table 1 also shows the relative importance of the quasi- 
static geometric errors, compared to the two other sources 
for the small feed rates. As expected, the summed con- 



tributions of the contouring errors and the dynamic errors 
are increasing from 1.7% to 14.3% for programmed feed 
rate of 1 000 mm/min and 18 896 mm/min respectively. 



Table 1: Repartition of the mean value of the error sources along the 
experimental trajectory for the six tests. 



F 
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1.1 


86.9 


11.4 


0.0 


0.6 


1800 


1.7 


85.9 


11.3 


0.4 
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3 240 
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84.7 


11.1 


0.9 


0.9 


5 832 


3.9 


82.6 


10.8 


1.1 


1.5 


10 498 


6.5 


77.7 


10.2 


2.3 


3.3 


18 896 


9.5 


72.8 


9.5 


3.3 


4.8 
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Table 2: Maximum errors of the three Cartesian components of the 
geometrical errors contributions in the 3 directions of the machine 
frame. 

Contribution x y z 

Link errors |<5 ; | ma:E (/jm) 
Motion errors |5 m | mo;E {um) 
Thermal drift \5td\max {um) 



Table 3: Maximum errors of contouring errors and dynamic errors 
contributions in the 3 directions of the machine frame for the 6 tested 
programmed feed rates. 
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1000 


8.1 


7.5 


10.1 


1.1 


1.7 


3.1 
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7.3 


11.2 


10.3 


1.5 


2.9 


3.2 


3 240 


11.0 


13.9 


10.6 


1.2 


3.5 


5.4 


5 832 


11.2 


16.9 


10.5 


1.8 


4.1 


7.1 


10 498 


11.5 


17.0 


12.2 
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6.2 


8.9 


18 896 


11.5 


17.0 


12.5 


3.7 


7.6 


10.4 



5.2. Maximum errors 

The ISO standards on geometric specifications [18] gen- 
erally define tolerance zones in which all the machined 
surface must be located. Considering this, the most pe- 
nalising errors are the maximum errors. 

The maximum errors |<5| ma a; of a contribution 8 are 
defined by the maximum absolute value among the terms 
of each column x 8, y 8 and z 8 of 8 (eq.(22|). 
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(max(| x <$|),max(| y 5|),max(| z <5|)) 



(22) 



Table 2 gives the maximum errors for the quasi-static 
geometric errors along the experimental trajectory. The 
maximum drift errors given are the maximum observed 
during the seven tests. Also, the effect of link errors de- 
pends on the pose of the machine, and the maximum error 
can be calculated offline for any tool path, and can be 
higher than the one observed in this study for other points 
in the working volume of the machine. 

Table 3 shows the evolution of the contributions of the 
contouring and dynamic errors with the programmed feed 
rate. The dynamic errors are not negligible, particularly 
when the programmed feed rate is over 10 000 mm/min. 
Even for the first test at F = 1 000 mm/min, dynamic 
errors occurred, but as shown in Fig. [9] only locally at the 
sharp corner of the experimental trajectory. The higher 
the feed rate, the higher the dynamic errors magnitude. 
As the feed-rate increases Fig. [10] shows that the dynamic 
errors also increase. 

5.3. Contributions along the trajectory 

Fig. [9] and Fig. [10] are two examples of the calculated 
contributions of 8 C , 8 qs and 8a along the experimental tra- 
jectory, respectively for the first and the sixth test. The 



Table 4: Evolution of the root mean square values of the errors in 
each direction for different feed rates. 

F 8d,rms {urn) 

{mm/min) x y z 
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3.4 
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Fig. 11: Evolution of Sdr-ms function of the programmed feed rate 
F in a logarithmic plot. The dashed line are curve fitted logarithmic 
models. 



acceleration of each axis is plotted to show that the con- 
touring errors and dynamic errors can, up to a certain 
point, be related to the acceleration. 

The block transition extracted from the NC unit is also 
plotted. The 17 linearly interpolated segments of the tra- 
jectory, previously mentioned in |4.3[ are visible on those 
figures. The expected acceleration peaks, which are re- 
lated to the block transitions, are visible in Fig. 



10 Those 



acceleration peaks are responsible for the sudden appear- 
ance of dynamic loads, themselves related to dynamic er- 
rors. In the conditions of those experiments, it shows that 
the machine geometry is affected by accelerations mainly 
due to the linear interpolation format. The peak-to-peak 
magnitude of the resulting dynamic errors is higher than 
5 fim along the y direction. 

The corner visible near the middle of the trajectory 
in Fig. [8] requires high acceleration of several axes, visi- 
ble in Fig. [9] and Fig. 10 The high dynamic load gener- 
ated on purpose results in visible dynamic geometric er- 
rors. Therefore, Fig. [10] shows that the machine accuracy 
is not equally affected in all directions. For instance, high 
acceleration of the X-axis seems to have little effect on 
dynamic errors according to Fig. 10 
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5.4- Root mean square of the dynamic errors 

The root mean square value of the dynamic errors along 
the trajectory in each Cartesian direction are given in Ta- 
ble 4. Fig. 11 shows the evolution of Sd^ rms as a function of 



Fona logarithmic plot. The dashed lines have equations 



of the form given by eq. ( 23 ) 



Orm.s — K ' r 



N 



(23) 



The values of k and N are obtained by curve fitting for 
each error and each direction. The aim of the curve fitting 
is not to propose a general model for dynamic RMS errors 
evolution, but only to notice a particular result for the ex- 
perimental conditions of this work. As the k and N values 
depend on the experimental trajectory, those results can 
not lead to general conclusion about the performance of 
the machine in terms of dynamic stiffness, but it suggests 
a new performance indicator under dynamic load with a 
standard trajectory designed for this purpose. 

6. Conclusion 

The high-speed machining context requires high accel- 
erations from the machine. Under those dynamic loads, 
the machine structure may no longer follow a rigid body 
behaviour. The new method presented in this paper al- 
lows to measure the volumetric error at the tool tip, and 
to evaluate the contribution of the dynamic geometric er- 
rors. 

The experiments at high programmed feed rates have 
been made possible by the use of a non-contact measuring 
instrument: the CapBall. The CapBall was re-designed to 
increase its stiffness for more reliable measurement under 
high dynamic loads. 

The evaluation of the respective magnitude of the con- 
touring errors, the quasi-static geometric errors and the 
dynamic geometric errors, using the mean values of the 
norm or the maximum errors, quantifies the relative im- 
pact of the dynamic errors, which can reach nearly 5% of 
the total volumetric error, considering the mean value of 
the norm. 

Furthermore, the study of the root mean square values 
of the errors shows that in this particular case, the evo- 
lution of the RMS errors with respect to the programmed 
feed rate can be approximated by a logarithmic law. 

Finally, the main interest of this experimental work is 
to propose a method to evaluate dynamic errors directly at 
the tool tip. Even if the method can only be applied on five 
- or more - axis machine tool, due to the close kinematic 
chain principle, it can be a powerful mean to validate mod- 
els for the dynamic behaviour of the structure with in- situ 
measurement. It has been shown that the linear interpola- 
tion generates acceleration and dynamic errors peaks. The 
influence of the NC interpolation and controller command 
law on dynamic errors can also be investigated with this 
method. 
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